Effect of Spin Correlations on Multi-orbital Metal-Insulator Transitions and Suppression of Orbital 

Selective Mott Transitions 
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We present the influence of spin correlation on the metal-insulator transitions (MIT) in two-orbital Hubbard 
models by the Kotliar-Ruckenstein slave-boson approach. In the asymmetric half-filling situation, the two orbits 
simultaneously transit from conducting to insulating states with the increase of Coulomb correlation, accompa- 
nied by a paramagnetic (PM)-antiferromagnetic (AFM) transition. The orbital selective Mott transition found 
in the PM condition is completely suppressed over a wide correlation range, though it may exist in the systems 
away from half-filling. In the insulating state, the system crosses over from a partially-polarized spin-gapped 
phase in the intermediate correlation regime to an almost fully-polarized Mott insulating phase in the strong 
correlation regime. These results demonstrate that the spin modulation to the quasiparticle spectra brings much 
rich and more interesting MIT scenario in multi-orbital correlated systems. 
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PACS numbers: 71.30.+h,75.30.Kz,71.10.Hf 

Orbital selective Mott transition (OSMT), i.e., with the in- 
crease of Coulomb correlation, a two-orbital electronic sys- 
tem transits from metal to a particular phase with one orbit 
exhibiting metallic conduction and another orbit being insulat- 
ing, is an interesting topic' . Theoretically it has been known 
that such an unusual phase is robust in asymmetric two-orbital 
Hubbard models^^. The OSMT phase can be driven by the 
asymmetric factors of two orbits, for example, the level split- 
ting between two orbits and the different bandwidths or band 
degeneration. There exist great debates on whether some typ- 
ical compounds with metal-insulator transitions (MIT), such 
as Sr2-.iCavRu04 and Va-.vCrvOa, are the prototypes of the 
OSMT phase^"^. This naturally arises a question: whether 
does the OSMT phase really exist in realistic compounds, es- 
pecially in the magnetic ordered systems? Generally speak- 
ing, the Mott-Hubbard MIT are accompanied with the change 
of magnetic structures. Such a change would crucially modu- 
late the quasiparticle states, hence significantly affect the MIT. 
The MIT and OSMT phase have been studied widely based 
on the dynamical mean-field theory (DMFT)** and Kotliar- 
Ruckenstein slave boson (KRSB) methods^, both of which are 
effective many-body approaches to treat the electronic inter- 
action over a wide correlation strength. However, limited by 
the huge computation time in the DMFT and KRSB methods, 
most of the studies focus on the paramagnetic (PM) phases. 
How the spin degree of freedom affects the MIT phase dia- 
gram and whether the OSMT phase exists in magnetic phase 
diagram are still not clear. 

Hesagawa'" and Sigrist et al.^^ first extended the single- 
orbital KRSB theory to the two-orbital situation to study the 
PM Mott transitions in two-orbital correlated electron sys- 



tems. To explore the influence of the spin degree of freedom 
on the Mott transitions in multi-orbital Hubbard models, we 
apply the extended KRSB method** on asymmetric two-orbital 
Hubbard models for various electron fillings. With the help of 
a new numerical ansatz to the multi-orbital KRSB solution 
we developed recently--, which overcomes the convergency 
problem of many parameters in minimizing the groundstate 
energy, we could treat arbitrary interaction strengths, finite 
interorbital hoppings or hybridizations, and various magnetic 
configurations in multi-orbital coiTelated electron systems. In 
this Letter, starting with an asymmetric two-orbital Hubbard 
model, we present various magnetic phase diagrams of the 
MIT, and discuss the physical scenario of the MIT in the pres- 
ence of spin degree of freedom. It is interestingly found that 
the OSMT phase completely vanishes at half filling. A Slater- 
type spin-gapped insulator and an orbital polarized band insu- 
lator are also found in magnetic phase diagrams, in addition 
to the conventional PM and Neel antiferromagnetic (AFM) 
metallic phases. Away from half filling, an AFM OSMT phase 
may be stable. These results show that the magnetic MIT pic- 
ture in multi-orbital correlated systems is more interesting and 
much richer than that without spin. 

Our starting point is an asymmetric two-orbital Hubbard 
model^ '*' applicable for various correlated electron systems 
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on a square lattice with only the nearest-neighbor hopping in- 
tegrals. Where tap, U iU') and Jh represent the hopping inte- 
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FIG. 1 : Magnetic phase diagram of the two-orbital Hubbard model 
in the n-U plane. and Ud denote the PM metal- AFM metal and 
the AFM metal- AFM insulator boundaries. Ud separates different 
orbital polarization regions. Theoretical parameters: D2/D1 = 0.5, 
Jh = 0.2(7 and A = 0. 



gral for the orbits a and /3, the intraband (inter-band) Coulomb 
repulsion and Hund's rule coupling, respectively. Throughout 
this paper we set U'- U - 2Jh, and take the energy band- 
width of the first orbit, 2Di = 2zfii - 2zt, be the energy unit. 
To explore the effect of spin in the multiorbital MIT, we in- 
troduce a few of auxiliary boson field operators representing 
the possibility of various electron occupations, such as e, p„o-, 
da-„a-'^, baj3, tacr and q, which denote various possibilities of 
none, single, double, triplicate and quadruplicate electron oc- 
cupations. Projecting the original fermion operators into these 
boson field and fermion field operators, one can obtain an ef- 
fective Hamiltonian, from which one can get the groundstate 
energy in the saddle point approximation. We apply the gen- 
eralized Lagrange multiplier method to enforce the fermion 
number constraint condition''^. In this framework the interor- 
bital hoppings and crystal field splitting can be treated in the 
same foot. 

We adopt an optimizing method to get the minimized 
groundstate energy for four different magnetic configurations, 
including the PM or nonmagnetic, ferromagnetic (FM) and 
Neel AFM states. In the optimizing process the boson field 
normalization set a boundary for the mean value of each bo- 
son field. As shown in our previous work'^, the optimizing 
problem of boundary constrained condition is still difficult. 
In numerically searching the global minima of the ground- 
state energy, we combine the pattern search method, the gra- 
dient method and the Rosenbrock method. As the optimizing 
point is approaching to the boundary, one must move one step 
inward the high-dimensional ellipsoid and the equipotential 
plane. Since the first axis of the new local orthogonal coor- 
dinate system of Rosenbrock method direct to the negative 
gradient direction, such an algorithm can be easily achieved. 
Our ansatz obtains a lot of desirable results for various PM 
systems with single, two and three orbitals'^ 

FiglUis magnetic phase diagrams of the asymmetric two- 



Q 0.6 

s 



(a) 0.0 





-■- (/, A=0.l,DyD=0.3 ■ 




-.- A=0.1,flj/£l =0.5 - 




Spin gap AFI 














PM Metal 





0.00 0.05 0.10 0.15 

JJU 



0.20 0.25 



Q 0.6 
<N 

s 



(b)o 



J -02V, DJD-0.1 
(/, J..-0.2U, DID -0.5 



AFI 



PM Metal 



0.00 0.04 0.08 0.12 0-16 0.20 

AI2D, 



FIG. 2: Magnetic phase diagram of half-filled two-orbital Hubbard 
model in the U-Jh plane with t^l2D\ =0.1 (a) and in the U-A plane 
with J„=0.2U (b) for the bandwidth ratio D2/D1 =0.3 and 0.5. 01 and 
AFI denote the orbital insulator and the AFM insulator, respectively. 



orbital Hubbard model with the bandwidth ratio D2lD\ - 0.5, 
the Hund's rule coupling Jh = 0.2U and the crystal field split- 
ting A = £2 - £1 =0. We only plot the phase boundaries for 
the electron filling of n < 2, the magnetic phase diagram for 
n > 2 can be mapped by considering the particle-hole sym- 
metry character in the present square lattice with the nearest- 
neighbor hopping. On the whole, we find four different sta- 
ble magnetic phases in Fig.l, including a PM metal, an AFM 
metal, an AFM OSMT phase and an AFM insulator The 
boundaries between these different phases are Uci and Uc3. 
Out of our expectation, the PM OSMT phase is not found for 
all of the electron filling and correlation strength. We address 
this magnetic phase diagram in length in the following. 

Firstly, we find that in the half filling (n=2) and A=0, with 
the increase of the electron correlation the system undergoes 
transitions from a PM metal (only for U=0) to an AFM metal 
at the the PM metal- AFM metal transition boundary Uci, and 
then to an AFM insulator at the AFM MIT boundary Uc3, 
as seen in Fig.l. It is obviously that the critical line Ud 
is considerable smaller than that without the spin degree of 
freedom'"*. Contrary to the early and recent results without 
the spin coiTelation''*- , the OSMT phase is completely sup- 
pressed in the present situation. Since the spin degree of free- 
dom is taken into account, the present AFM insulator is a spin 
gapped state in the intermediate correlation region, or a AFM 
Mott state in the large U region. Approximately, the spin 
gapped insulating region at half filling covers the PM OSMT 
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phase region obtained by the DMFT approach^'*. This shows 
that when U increases, the narrow orbit of the system becomes 
of AFM insulating. Such an AFM correlation also modifies 
the electronic properties in the wide band, and leads the spec- 
tra of the wide-band electron to open a gap. Hence the MIT 
simultaneously occurs in the two-orbital system. As a result, 
the PM OSMT phase is completely suppressed at half filling. 

Away from half filling, the magnetic phase diagram be- 
comes of more interesting. As seen in Fig.l, though the PM 
OSMT is completely suppressed, the AFM OSMT phase can 
survive over a wide region, in which the half-filled narrow 
band becomes of AFM and insulating, while the wide band 
away from half filling remains conducting, though modulated 
by the AFM correlation. We find that with the increase of the 
Coulomb correlation, the system undergoes a PM metal- AFM 
metal transition and an AFM metal-AFM OSMT transition. 
These transitions are also separated by the critical line Ud and 
Uc3, respectively. Even in the PM phase, the system may show 
different orbital polarization separated by the critical line Vd, 
as indicated in Fig.l. When U < Ud, the particles are mainly 
distributed in the wide band and the reverse particle distribu- 
tions achieve when U > Ud, since the Coulomb interaction 
drives a charge redistribution with increasing t/Ji^i. We no- 
tice that the critical line Ud in Fig. [T]is almost independent 
of the magnetic transition, suggesting that the charge redis- 
tribution in different orbits is directly driven by the electron 
correlation, rather than the spin exchange splitting. 

Though the PM OSMP is suppressed in the presence of 
magnetic correlation at half filling and A = 0, it is unknown 
whether it exists over wide crystal field splitting A and the 
Hund's rule coupling. In what follows, we focus on these sit- 
uations in the half-filled Hubbard model. The magnetic phase 
diagrams for the half-filled asymmetric Hubbard models in the 
U - Jh and U - A planes are plotted in Fig|2] It is still find 
that the OSMT phase is absent in Fig. 2a and Fig. 2b. Also, as 
seen in these two figures, the AFM metallic phase is unstable 
for finite crystal field splitting, it will become of a PM metal- 
lic one. Thus increasing the crystal field splitting may drive a 
AFM-PM transition. 

Fig|2^ shows the U - Jh phase diagram for different band- 
width ratios. In the PM metallic situation below the phase 
boundary Uc3, the electrons in two orbits are itinerant and 
equivalently occupy two orbits. The two orbits are singly oc- 
cupied by two electrons and become insulating and the PM- 
AFM MIT occurs as U > Ud ■ In the presence of spin degree 
of freedom, the system enters the spin-gapped AFM insulat- 
ing phase in the intermediate region in Fig. 2a. Due to the spin 
modulation, the PM-AFM MIT boundary Ud is considerable 
smaller than that without spin. It is interesting that when the 
two orbits are split by a crystal field, a new insulating phase, 
named orbital insulator, appears in the small Jh region. In the 
orbital insulating phase, the two electrons fully occupy the 
lower narrow orbits and the system is nonmagnetic and enters 
into an insulating phase. Such an interesting phase has also 
been obtained recently by Yu and Si'^. Due to the neglect of 
the spin effect in Yu and Sis work, their Uc/2Di is larger than 
unity, considerably larger than our result, about 0.68, in the 
system with Jh = 0, D2/D1 = 0.3 and A = 0.1. With the 



further increase of the electronic correlation, one may expect 
that the system crosses overs from a spin-gapped insulator to 
a conventional Mott insulator, as we will show in Fig. 3 later. 
These unusual results show that the spin correlation has pro- 
found influence on the quantum phases and the boarder of the 
MIT in correlated electron systems. 

On the other hand, the magnetic phase diagram as the func- 
tion of U and A is shown in Fig|2]3. One can see that the 
OSMT phase is also excluded. From the U-A phase diagram in 
Fig|2j3, we see that the critical line Ud of the PM-AFM MIT 
decreases with the decrease of the bandwidth ratios D2/D1. 
At the half filling, large Coulomb correlation and Hund's rule 
coupling, which contribute the spin exchange splitting, are 
needed to overcome the energy increase arising from the crys- 
tal field splitting. So the MIT boundary Ud lifts with the in- 
crease of A. It is interesting that the spin-gapped AFM insu- 
lator in Fig.l and Fig. 2 is not observed in the previous litera- 
ture. In this phase the two orbits are not fully spin polarized, 
this allows spins to form density wave order Considering the 
fact that this insulator is stable in the intermediate correlated 
regime, we attribute it to the Slater-type insulator'^ with a par- 
tial spin polarization. 
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FIG. 3: Correlation dependence of occupation probabilities (a), 
band renormalization factor (b), particle number (c) and density of 
states (DOS) (d) of the half-filling two-orbital Hubbard model with 
O2/O1 = 0.5, Jh = 0-2U and A/2Z), = 0.1. Inset in (a) displays the 
full data of the boson occupation probability. 

Figl3]shows the correlation dependence of occupation prob- 
abilities, band renormalization factor, particle number and 
DOS for finite Hund's rule coupling in the half-filled two- 
orbital Hubbard model. All these figures clearly show that 
the MIT occurs at the critical point Ud, accompanied by a 
magnetic transition. Fig|3^ displays the U dependence of var- 
ious occupation probabilities. It shows that a first-order- type 
MIT transition occurs at Ud - 0.42: when U < Ud, various 
charge and spin configurations appear in approximately equal 
weight, suggesting that the PM metallic ground state is stable. 
When U > Ud, the spin-triplet dg-cr occupation is dominant, 
the system enters the magnetically ordered insulating phase. 
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As seen from Figl3^, only the dn occupancy with each orbit 
occupied by one electron increases with the lift of U and Jff, 
the other occupancies decrease immediately when U > U^. 
After U >- 2D\, all of the occupancy probabilities almost 
approach constants. This indicates that the system enters a 
conventional Mott insulating phase. 

FigOh shows the dependence of the band renormalization 
factor Z„ of the quasiparticles in orbit-ff on the Coulomb cor- 
relation U, which significantly differs from the behavior of Z 
without considering the spin correlation effect. In accordance 
with Fig|3^, we find that with the increase of U, the renormal- 
ization factor decreases in the PM phase since the electronic 
correlation narrows the bandwidthes. Accompanied with the 
transition from the PM metallic phase to the AFM insulating 
one, an insulating spin gap appears. Out of our expectation, it 
abruptly increases at the MIT critical point Uc3/2D-i = 0.42, 
and then gradually increases abnormally. With the increase 
of U/2Di to about 1.0, both Zi and Z2 approach saturations. 
The system may cross over from a Slater-type insulator to a 
conventional Mott-Heisenberg insulator, though both of the 
band renormalization factors do not show such a crossover 
Our results are distinctly different from the PM results by 
Hesagawa'" and Sigrist et a/. " , where there exists a stable and 
robust OSMT phase in the intermediate correlation regime. 

The U dependence of the particle number of each orbit per 
spin channel is shown in Fig[3};. It is clearly seen that with the 
increase of the electronic correlation, the particle numbers in 
the two orbits differ from each other in the PM region and ex- 
perience a critical change at Uc3, and gradually approach con- 
stant occupations as U jlDi approaches unity, in accordance 
with Fig.3a and Fig.3b. Obviously the PM-AFM MIT at t/^ 
is first order Although the spins are not polarized in the PM 
metallic phase, the sublattice spins are strongly polarized in 
the insulating phase ({/ > UcT,) and the sublattice spin polar- 
ization almost saturates as U/2D[ 1. Asymmetric band- 
widths and crystal field splitting lead the spin polarizations of 
two orbits be different. The more interesting is that with the 
transition from the PM to AFM phases when U > Ud, the 
increase of the Coulomb interaction depresses the orbital po- 
larization, which originates from the increase of the Hund's 
rule coupling. We find that the finite crystal field splitting is 
the key reason for the first-order MIT. When the crystal field 
splitting equals to zero, the transition from the PM to magnet- 



ically ordered phases is the second order, in agreement with 
Hasegawa's results'*^. 

The spin-dependent densities of states (DOS) for U/2Di = 
0, 0.37, 0.45 and 0.5 are plotted in Fig|3}J. It demonstrates how 
the two-orbital DOS evolve with the electronic correlation. 
In the PM metalUc phase, both bandwidths of the two orbits 
shrink with increasing U. In the present square lattice, one 
finds two cusps in the spin-dependent DOS of the two orbits 
when the crystal field splitting is switched on. The MIT en- 
ergy gap, which becomes large with the increase of U, opens 
at Uc3 - 0.42 in the system with D2/Di=0.5, Jh=0.2U and 
A/2Di=0.1. The well-known quasiparticle peaks in Fermi 
surface and the three-peaks structure found in the DMFT ap- 
proach are not observed in the spin-dependent multi-orbital 
KRSB theory, which is attributed to the neglect of dynamical 
transition process between Hubbard subbands in the present 
mean-field approach. 

From the preceding studies on the metal-insulator transi- 
tions and the orbital selective Mott phase in the two-orbital 
Hubbard model with spin degree of freedom in the slave bo- 
son mean field theory, we find that in the presence of spin cor- 
relation, the scenario of the metal-insulator transitions may 
be distinctly different from that in the paramagnetic situa- 
tions. In the half-filled asymmetric two-orbital Hubbard mod- 
els, the spin modulation arising from the insulating narrow 
band opens a gap in the wide band, leading the paramagnetic 
orbital selective Mott phase to be suppressed over a wide cor- 
relation range; and the orbital insulating phase is also found 
to stable in the small Jh region. In the presence of spin cor- 
relations and finite crystal field splitting, the metal-insulator 
transitions, accompanied by magnetic transitions, are of first 
order Away from half-filling an antiferromagnetic orbital se- 
lective Mott phase may survive and the critical point of U£-3 
is smaller than the PM results because the symmetry of the 
local interaction is lowered by the spin polarization. These 
results demonstrate that the spin correlations and the crystal 
field splitting play very important roles for the metal-insulator 
transitions. 
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